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Abstract 

A class of singular 3I?-velocity vector fields is constructed which 
satisfy the incompressible 3D-Euler equation. It is shown that such 
a solution scheme does not exist in dimension 2. The solutions con- 
structed are bounded and smooth up to finite time where they become 
singular. Although the solution is smooth and bounded there seems to 
be no bound in of the velocity field. 
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1 Statement of results and idea of construction of 
singular solutions 

The main results of this paper concern singular solutions of the n-D Euler 
equation for n > 3, i.e., the equation system 

dviisr^n dvj _ dp 
at "I" Z^i=l dxj ~ dxi ' 

EILi = 0, (1) 

Vi{Q,x) = hi{x). 

For the main part of this paper we consider evolutions on the whole domain 
of M"". Analogous results for the n-torus T" with respect to dual Sobolev 
spaces may be derived, but we did not consider the details so far. It is 
desriable to have a construction for data of physical interest which are C°° 
and such that all multivariate derivatives are of polynomial decay of order 
m > 2, i.e., for multiindices a = (ai, • • • , a„) with > and \a\ = X^ILi 
we have 
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where -D° = |^ denotes the multivariate derivative of order a. In the 
present paper we construct a bounded smooth vector field where the first 
two components are of finite energy. It may be interesting to consider similar 
constructions for derivatives of the value functions and possibly with weaker 
singularities in t- although the connection of data and solution may not be 
that simple as in this paper implying that the construction has to be more 
involved. If a function / is of polynomial decay of any order m S N then we 
simply say that h is of polynomial decay. Polynomial decay and smoothness 
implies that the data are in the intersection PlseR-ff'* of standard Sobolev 
spaces i/* = (M") of order s G M, and hence polynomial decay of the 
Fourier transforms. On the n-torus dual Sobolev spaces may be used. In 
this case multiindexed vectors 

(the superscript T meaning 'transposed') may be considered with complex 
constants which decay fast enough as |a| t oo corresponds to a function 

n G C°° iXi) , (4) 

where 



) := ^ -Ua exp { ^^^"^ j . (5) 



We say that the infinite vector u'^ is in the dual Sobolev space of order 
s G M, i.e., 

G (Z") (6) 

in symbols, if the corresponding function u defined in Q is in the Sobolev 
space (Tp) of order s G M. We may also define the dual space directly 
via 

G (Z") ^ \uj^ {af' < oo, (7) 

where 

{a):={l + \a\Y\ (8) 

The two definitions are clearly equivalent. 

Our first observation reduces the search for a class of singular solutions 
to the search for a class of nonsingular global solutions of certain nonlinear 
partial integro-differential solutions. We have 

Theorem 1.1. // n = 3 then for all data /i,/2 G if^ n and /a G 
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which satisfy the nonlinear partial integro- differential equation 

- fl + + /2/l,2 - /k3 K^X - y) (/f i + + + /2,2)' + 



/l,2/2,l + /l,3^3,l(/l,l> /2,2) + /2,3-f3,2(/l,l, /2,2)j iy)dyjf2,3 

= ( - /2 + /l/2,l + /2/2,2 - / K,2{X - y) (/2 , + /f ^ + + /s.s)' + 

/l,2/2,l + /l,3^3,l(/l,l'/2,2) +/2,3-f3,2(/l,l,/2,2))(?/)dy)/l,3, 



(9) 

i/ien there exists a singular solution to the incompressible Euler equation 
which becomes singular at finite time to > 0, and is smooth for t < t^. Here, 

fsix) = /3(/i,i,/2,2) := r - /2,2) {xi,X2,y)dy (10) 



A similar result holds for n > 3. 

Remark 1.2. A similar result may hold on the n-torus. 

A second observation is that solutions to ([9]) exist, and in a rather generic 
sense, i.e., in a sense which makes it impossible to single out the data set of 
the constructed set as 'exceptional'. We have 

Theorem 1.3. For each fi £ H'^CiC'^ which is of polynomial decay of order 
2 there exist f2 S H of polynomial decay of order 2, and f^ G such 
that the vector field {fi, f2, fi)'^ is a global solution of equation More 
general, for or each f\ € n C™ which is of polynomial decay of order m 
there exist f2 G -fT™ H and fs € of polynomial decay of order m such 
that the vector field (/i, /2, f^)"^ is a global solution of equation 

Next let us consider some background for the results above. The Euler 
equation is an equation without diffusion. If we add an diffusion operator 
—u/S. with a constant viscosity > we get the incompressible Navier-Stokes 
equation 



n=i ^ = 0, (11) 

Vi{0,x) = hi{x). 

If we cancel the pressure p and incompressibility we get the multivariate 
Burgers equation. 

(12) 

Vi{Q,x) = hi{x). 
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Here the gi G Ds^rH^ are additional source terms which may occur also in 
the formulation of the incompressible Navier-Stokes equation. For u > 
global existence for the multivariate Burgers equation is known. An ele- 
mentary proof may be found in [3] and [1]. For the multivariate Burgers 
equation it is more or less obvious that global existence is not preserved 
as the equation degenerates, i.e., as J, 0. In this case, the function 
Vi-.W^ ^R, l<i<n defined by 

v,it,x):=^, (13) 



along with 

dfijx) 
dxj 



i=i 

for all 1 < i < n and all x G M" solves the inviscid equation 

dvi dvi 

for 1 < i < n, and with initial data 

Vi{0,x) = hi{x) =: -fi{x). (16) 

Obviously, this velocity field becomes singular at t = 1 as long as /j 7^ 
for some 1 < i < n. This is a phenomenon which we can observed in all 
dimensions n > 1. We may say that viscosity = is a bifurcation point 
for existence of global behavior. Indeed for < the equation is even 
ill-posed in most situations. This situation becomes more complicated and 
more interesting if we consider the analogous relationship between the in- 
compressible Euler equation and the incompressible Navier-Stokes equation. 
For in this case we know that that the Euler equation has global solutions in 
dimension n = 2 and we maintained that the incompressible Navier-Stokes 
equation has global solutions in different situations (cf. 0], [2]j and [1]). The 
difference of the Euler equation and a degenerate Burgers equation is the 
existence of pressure and incompressibility, and the relation of both which 
yields elimination of the pressure by the Leray projection. If the divergence 
of the vector field v = (t;!, • • • , Vn)'^ of the form (fT3]l is zero, i.e., if 

n r) n 
i=l * 1=1 

then this incompressibility together with the relation ()14p implies that the 
divergence of the vector field f = (/i, • • • , /n)^ is zero, and whence 

- EILi M^) + EILi E •=! (/.(^)/m(^)), 

(18) 
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However the Leray projection form of the incompressible Euler equation is 
just 

(19) 

Vi{0,x) = -fi{x). 

Here, denotes the fundamental solution of the Poisson equation in di- 
mension n and K^^i denotes its partial derivative with respect to the variable 
Xi. This leads is to the conclusion 

Proposition 1.4. The function Vi : M" — )• M, 1 < i < n defined for all 
1 < i < n by 

v,{t,x):=-^^ (20) 

satisfies the incompressible Euler equation with initial data hi = —fi : M" — )■ 
M if the conditions 

- m + f2 ^^(^^^^ = 0, 1 < i < n, (21) 



j=i - - J 



and 



and 



^f,4x)fij{x) = 0, (22) 

n 

divf = Y, ki = (23) 

i=l 



are satisfied. 



If fi satisfy periodic boundary conditions, then an analogous statement 
is true on the n-torus of course. For a torus of size / the relations (|14p may 
be written in the basis {exp (^^^^) and where kx = ^kiXi along 

with X = (xi, • • • , Xn)- This leads to the relation 

n 

/ia = - ^ ^ fj(^a-f3)Pjfil3 (24) 



for the modes fia for 1 < i < n and the multiindices a, P € TP" Next to 
the two observations stated at the beginning of this section we make a third 
observation which we would expect from the fact that global solutions to 
the Euler equation exist in dimension 2. We have 

Theorem 1.5. In case of dimension n = 2 the only classical solution of 
form [2U\) is 

fi = (25) 

for all 1 < i < n. This is compatible with the well-known fact that classical 
global solutions for 2-D incompressible Euler equations exist. 
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In the next section we prove this theorem. Then in section 3 we prove 
theorem 11.11 and in section 4 we prove theorem II. 3[ In the last section we 
prove the existence of global solutions to a Cauchy problem which is related 
to the partial integro-differential equation of theorem II. 1[ We may summa- 
rize the considerations of this and the mentioned related papers saying that 
viscosity i/ = is a bifurcation point of global solution behavior in case of 
dimension n > 3. 

2 Proof of theorem 11.51 

We prove the theorem on the n-torus T" for n = 2. For dimension n = 2 
divergence becomes 

ViA{t,x) = -V2,2{t,x) (26) 

for all (t, x), hence 

= -h,2{x). (27) 
Using this incompressibility, and differentiating the relations 

Vi + f 1^1,1 + V2Vi^2 = 

(28) 

V2 + f 1^2,1 + V2V2,2 = 

with respect to xi and X2 respectively, summing up, and multiplying by 
(1 — t) leads to 

+ 2/1,2/2,1 + /|,2 = 0. (29) 
Along with divergence the latter equation leads to 

/l,2/2,l + /2,2 = /l,l + /l, 2/2,1 = — /l,l/2,2 + /l, 2/2,1 = 0. (30) 

The last equation of (|3Up may be rewritten in the form 

det ( {^'1 \ = 0. (31) 

V /1,2 /2,2 / 

Hence, the vectors (/i,i, /i,2)"^ and (/2,i,/2,2)^ are linearly dependent, i.e., 
there is a A G M such that 

/i,i = A/2,1 and /i,2 = A/2,2- (32) 

The latter relation may be considered together with the incompressibility 
relation 

/l,l = -/2,2. (33) 

Consider the ansatz 

= E -^^^P (34) 
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for 1 < i < 2, and a = (ai,a2) G and ax := a\X\ + a2X2- Note that 
psp is equivalent to the statement that for all a G 



/a"! = -/aa2- (35) 

This implies that A = — 1 in (j32p . Now from (j28p we have 

-/l + + /2/l,2 = 

(36) 

-h + /l/2,l + /2/2,2 = 0, 

and this together with the information A = — 1 in ()32p and incompressibility 
leads to 



(37) 



-/l + + /2/l,2 = -/l + - /2/2,2 = 0, 

-/2 + /l/2,l + /2f2,2 = -/2 - /l/l,l + 72/2,2 = 0, 

Adding both equations we get 

h = -h- (38) 
This together with ()37p and vanishing divergence leads to 

-/l + - /2/2,2 = -/l + /l(/l,l + /2,2) = -/l = 0, 

(39) 

-h - + /2/2,2 = -h + h{h,l + /2,2) = -h = 0, 

and we are done. 

3 Proof of theorem 11.11 

Let us first explain why the theorem is stated with finite (maybe smaller) 
io > instead of to = 1 as may be expected from reading the considerations 
made so far. Indeed, in the following we define 

^.(t,x):=^ (40) 

and derive an equation as in the statement of theorem 11.11 above. However 
this imposes the condition 

/l,l + /2,2>-l (41) 

which is satisfied for sufficiently small i/^-norm of /i and /2. However 
reading the argument below with the vector field 

fix) 

Vi{t,x) := j^^, 1 < i < n, (42) 
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leads to the concusion that the restriction in ()4ip is no essential restriction. If 
we choose such a vector field, then the argument below imposes the condition 

+ /2,2 > -c, (43) 

and choosing c > large enough this is no essential restriction: given some 
/i € H'^nC'^ we shall show that there are /2 G -ff^nC^ and /s G such that 
we have a divergence free vector field, and such that the pair (/i, /2) solves 
the partial integro-differential equation of theorem 11.11 We may then choose 
c > as an upper bound for |/i,i| + |/2,2| and apply the following argument 
in order to show that for this data we have a solution which becomes singular 
for some to < K The reader should think of small i/^-data first such that 
the condition in ()4ip is satisfied and then consider the preceding remark in 
order to convince himself that the condition (|4ip is a convenient assumption 
which can be easily dropped. We prove the theorem for and make some 
remarks that it applies also to the case where the domain is the n-torus for 
n = 3 as we go along. We construct a vector field (/i, /2, /a)"^ with functions 
fi : C°°, 1 < i < n, and which satisfies the following three conditions 

i) The data satisfy a multivariate Burgers equation 

-fl + + /2/l,2 + /3/l,3 = Gi 

-f2 + /l/2,l + /2/2,2 + /3/2,3 = G2 (44) 

— /3 + fih,i + 72/3,2 + 73/3,3 = G3 
where Gi G are some source term functions chosen below. 



ii) The divergence of the vector field is zero, i.e., 

+ /2,2 + /3,3 = (45) 

iii) The Leray projection term of the vector field (/i, /2, fs)^ satisfies 

+ /l,2 + /|,3 + /l,2/2,l + /l,3/3,l + /2,3/3,2 = 9 (46) 

for some function g £ C'^. In case of the domain M" the condition 
takes the form 



Gi{x) = j K4x-y)g{y)dy, (47) 

where K is the Poisson kernel in dimension = 3 and K^i denotes 
its partial derivative with respect to the variable Xj. Furthermore in 
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case of the n-torus g is such that the functions Gi above satisfy for 
1 < i < n 

— v=^n 2 ' 

2^1=1 "^i 

where 

, , v-^ / 27riax\ , , , 

f (2;) = l^9oc exp — - — (49) 

(In the following we consider the case / = 1 w.l.o.g..) 
It is clear then that the velocity field 

Vi(t,x)=^-^^, \ <i<n (50) 

satisfies the incompressible Navier Stokes equation in its Leray projection 
form with initial data {fi{x), f2{x), /^{x))'^ on the domain [0,1) x M" (and 
[0,1) X T" respectively). Equivalently, it satisfies a multivariate Burgers 
equation 

Q. 

Vi,t + VlVi,l + V2Vi^2 + V3Vi^3 = J^Z^' ^^^^ 

for 1 < i < n by construction. We want to show that for /i G n 
there exists f2 G H and /s G such that conditions i),ii), and iii) 
are satisfied. First for any data /i,/2 G n we have decay to zero 
at infinity for derivatives up to second order, i.e., as |x| t 00. Hence from 
incompressibility (condition ii) above) we have 

{-fi,i-f2,2){xi,X2,y)dy (52) 

-00 

(we may invoke the assumption of some polynomial decay here but it is not 
necessary). Next we reduce condition i). The third equation of condition i) 
together with incompressibility /3^3 = — — /2,2 leads to 

-h + hhl + /2/3,2 + h (-/l,l - /2,2) = G3. (53) 

This means that 

h = 7 7 , (54) 

-J- - - 72,2 

where + /2,2 / 1 and 

/l/3,l - /2/3,2 = G3 (55) 

otherwise. For the moment we assume that 1 is not in the range of + /2,2 
such that the first alternative holds for all arguments x G T". Then condition 
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i) reduces to 



-fi + + /2/l,2 + \ '\-^fi,i^ffJ^'^ ) /i.3 — Gi 



-f2 + /l/2,l + 72/2,2 + ( '^'^ l-^fi i-ht '^ ) -^2,3 — G2 



(56) 



We have not defined Gi, G2, and G3 yet. They are defined via the condition 
iii). We are free to define g in this condition, i.e., we define 

9 ■= 7i,i + 71,2 + 71,3 + 71,2/2,1 + 71,3/3,1 + 72,3/3,2 

(57) 

= 7l,l + /2,2 + (/l,l + /2,2) + /l, 2/2,1 + /l,3/3,l + /2,3/3,2 

The system (j56|) is linear in f^^i and /3^2- We get 

(1 - /l,l - /2,2) i-fl + /l/l,l + /2/l,2 - Gi) = (-G3 + /l/3,l + /2/3,2) /l,3 
(1 - /l,l — /2,2) ( — /2 + /l/2,l + /2/2,2 — G2) = ( — Gs + /l/3,l + f2h,2) /2,3- 

(58) 

The first equation multiplied by /2,3 has to equal the second equation mul- 
tiphed by /i,3 

(1 - /l,l - /2,2) [ (-/l + /l/l,l + /2/l,2 - Gi) /2,3 

(59) 

- (-/2 + /l/2,l + /2/2,2 - G2) /l,3] = 0, 

and according to our temporary assumption that 1 is not in the range of 
/i,i + /2,2 this means that 

( — /l + /l/l,l + /2/l,2 - Gi) /2,3 = (-/2 + /l/2,l + /2/2,2 — G2) /l,3- 

(60) 

In the latter equation Gi, 1 < i < 2 are functionals of /i and f2- 
Gi{x) = J K^i{x - y)g{y)dy 

= J^3K^i{x-y)x 

/l,l + /2,2 + (/l,l + /2,2)' + /l,2/2,l + /l,3/3,l(/l,l, /2,2) + /2,3/3,2 (/l,l , /2,2) ) {y)dy 

(61) 
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Hence we get the integro-differential equation 

' fl + + /2/l,2 - K^iix - y) (/f i + + + /2,2)' + 

/l,2/2,l + fl,3h,l{fl,l, f2,2) + /2,3-^3,2(/l,l,/2,2))(?/)dy)/2,3 
= ( - /2 + /l/2,l + /2/2,2 - / i^,2(x - 2/) (/f, 1 + /|,2 + + /2,2)' + 

/l,2/2,l + /l,3^3,l(/l,l'/2,2) + /2,3-^3,2(/l,l,/2,2))(?/)d?/)/l,3. 



(62) 

On the 3 torus a similar expression for Gi for / = 1 may be defined in 
terms of g by 

Eaezs ((27 rzai))ffae xp (27riaa;) 
and together with the representations 



"-^i — ^^n 2 ' V""^/ 



fiix)= f^exp{27riax), (64) 

and 

/2(a:) = ^ exp (27r«a3;) . (65) 

aeZ3 



4 Proof of theorem 11.3 

We rewrite ([9]) as 

{ — fl + + /2/l,2)/2,3 + (/2 — /l/2,l — /2/2,2)/l,3 

i^,l(x - y) (/f l + /1 2 + (/i,i + /2,2)' + 

/l,2/2,l + /l,34,l(/l,li /2,2) + /2,3-^3,2(/l,l,/2,2))(y)rfy)/2,3 (66) 
/ /i',2(x - y) (fl, + /|,2 + (/l,l + /2,2)' + 



/l,2/2,l + /l,34,l(/l,l> /2,2) + /2,3-^3,2(/l,l,/2,2)j(y)rfyj/l,3• 
In the following discussion we consider the case n = 3 and work on the whole 
domain M^. Similar considerations hold for the n-torus. Given /i S H'^CiC^ 
you may try to construct /2 as a limit /2 = limfci-oo /2 where for a start value 
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/2 = — /i you may try the iteration scheme 

((-/l + + /|7l,2)/|,3 + ift' - hfil - f2f2,2)fl,3 = 

/j,3((l - <Pi)K),,ix - y){{h,i? + {fl,? + + fl,f + 

/l,2/|,l + /l,3/3,l(/l,l, fi2) + f2,3hAfl,U fl2)) (y)%) fls (67) 

-[jK2{x- y) [fl, + {fl,f + + fl,) ' + 

/l,2/|,l + /l,3/3,l(/l,l, /|,2) + fysAfhU /|,2)) (2/)d2/)/l,3. 

The problem with this is that it is not easy to find an appropriate function 
space for the iteration, i.e., if you start an iteration step with /2 G H^, then 
it seems difficult to observe that f^^'^ G H"^. For this reason we consider the 
following Cauchy problem: 

/i,3 Ki{x - y) [fl, + (^^2)' + + <2)' + 

/l,2t''l + /l,3-^3,l(/l,l,vf2) + f2,3h,2{h,l^^%))iy)dy)v% 

< 

- ( ^3 K2{x - y) (/fi + [v^^^f + + v^^f + 

/l,2'f|,l + /l,3^3,l(/l,l'^2,2) + V%h,2{hA^'"%)){y)dy) fl,^^ 

v\Q,x) = 0. 

(68) 

Theorem 4.1. T/ie function : — )• M defined for all x G 6?/ 

fj^(x) := limlimf'^(t,x) (69) 

e4,0 t\oo 

is well defined, and such that G H C^. 

Proof. For e > let f ^ be the solution of the Cauchy problem (|68p which we 
construct in the following section. In terms of the fundamental solution of 

— - eAv^ = 0, (70) 
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which is (for t > s) 



exp ( -^-^ ) , (71) 



(2v/e7r(t- s))" V 4e(t-s)y' 



the solution f ^ of the Cauchy problem (j68|) has the pointwise representation 
(with s' = {t — s), therefore an additional minus sign) 

V^t, X) = ^3 ((-/l + + ^71,2)^^^3 + - fiyy - V'v%)h,^) 

X , i-, exp i-^^r^) dzds' 

+ ( ^3 K i(^ - 2/) (/I'l + {v^f + + <2)' + 

/l,2^;,\ + /l,3/3,l(/l,l. ^:2) + /2,3/3,2(/l,l, ^^:2)) (2/)%) ^3) ^) 

- /o ( /r3 i^,2(x - y) [fl, + {V^,,f + + V^^^f + 

/i,2^^2,i + h,?>hMl^^'"2,2) + v%hMl^^v^,2)){y)dyfl^^{s,z) 

X , i-^, exp f-^^) dzds. 

(72) 

The time derivative is the evaluation of the integrand on the right side of 
T2\f at time t. Note that this integrand evaluated at t has an upper bound 

(73) 



(2Ve7ft)" 

which goes to zero as 1 1 oo. It follows that 



lim— = 0. (74) 
ttoo at 



Furthermore, 



IS 



vl,:=\unv\t,.) (75) 

tfoo 

well defined, and we have G n according to the result on the 
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Cauchy problem of the next section. Moreover solves 

- ( /„3 K2{x - y) [fl, + iv%) ' + + v^^,) ' + 
fi,2v%i + fi,3h,i{f 1,1,^2^2) + '"%hAfhi^'"%)) {y)dy^f 1,3- 

(76) 

This holds for all e > where |v'^|/f2 n < C for some finite C > which 
is independent of e > (cf. next section) . Hence the limit /2 = lim^j^o 
exists and solves the equation ([76]) for e | 0. □ 



5 Proof of global existence for the Cauchy prob- 
lem 

We prove global existence via an iteration. There are several possibilities 
which part of the iteration of the linearization is put into a source term and 
which part is taken into the definition of the construction of the fundamental 
solution at each iteration step. In [T] and in [3] and [1] we used the adjoint 
of fundamental solutions in order to estimate the local growth of such kind 
of iterations in a H^-norm. In this paper we simplify the argument a little 
and avoid the adjoint and some estimates, i.e, we represent iterative approx- 
imations directly by convolutions with a Gaussian. An additional difficulty 
(compared to the schemes of the incompressible Navier Stokes equation) are 
the integral terms involving the operators I^^i and 13^2 • It seems that a 
contraction property for the functional increments of our scheme cannot be 
proved directly. However, these terms just ensure incompressibility of the 
vector field. Therefore we replace these integrals terms first by a param- 
eter function w{t, .) G and then solve the Cauchy problem for a fixed 
parameter w. Then we set up a second integral equation corresponding to 
incompressibility, and the solution to the latter equation is the solution leads 
to the solution of the Cauchy problem. We start with v"^'^ = — /i, and given 
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ye,w,k define v''''^^^^^ by the iteration 

V^'-'>'+\t,x) = ^3 (((-/l + + V^'-''fl,2)v:t' 

fi,2V^f'^ + f 1,3^^1 + v''^^'"''^w^2^ {y)dy^ ■^'3"'''') (s, 2)^, (t - s, X - z)dzds 
- So ( /m3 K2 {x - y) [fl + {v:f'') " + + <f ' + 
fh^v^'J''' + /i,3W^,i + I'^r'^'^.s) {y)dyf 1,3^ (s, z)Ge(t - s, X - 

(77) 

This is the solution scheme for a Cauchy problem parameterized by w{t, .) £ 
C|, where the corresponding naive iteration scheme would be one without 

diffusion. Note that Ge(t — s,x — y) = — , ^ exp ( — wT^\ | is the 

fundamental solution of 

^ - eAp = 0. (78) 
For u : [0, 00) x — ^ M we define 

j-uj^P'*^ := sup exp (-Ci) |u(i, .)|^2, (79) 
te[o,oo) 

where C > is determined below in order to obtain a global contraction. 
Next we prove a contraction property for the parameterized Cauchy problem, 
i.e., we prove for given w{t, .) G 

|j^.,^,fc+ijexp,c ^ l|^^,,^,fc|exp,c_ ^gQ^ 

In the following we denote the space associated with the norm (|79p with a 
specific C > by -ffgxp c- ^^xt for the functional increments (5ti'^'"''^'+i _ 
^e,fc+i _ ye,w,k have the Cauchy problem 
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dt 



eAdv' 



(81) 



^ J^;^.«',fc+i(0,2;) =0, 



where 



5^'+' = /M3i^,i(x-y)(/2i + 



t,w,k 



e,w,k 



+ 



e,w,k 



-{j^sK2ix-y){fl + 

fl,2v\'"''' + fl,3W,l + V%'"'^W^2 ) (y)c?y ) /l,3- 



(82) 



Note that the source term 5'^'*^+! for the original Cauchy problem is obtained 
by replacing w^2 by /3,2(/i,ii '^^2'') ^^id wi by /3,i(/i,i) '^^^2'') i'^ S'^''^"^- We shall 
construct the functions v'^''^'''{t, .) G H'^HC'^ uniformly with respect to t > 0. 
The limit of the iteration is then in C C° uniformly in t, and this leads to 
representations of the value function v'' which imply more regularity. First 
note that we can represent the solution Sv'''^''''^^ of the Cauchy problem ()8ip 
in the form 



,w,k+l 



/o /r3 ((-/i + hf 1,1 + ^^•'"'Vi,2)5<r'' + 6v^'-'''h,2v:t' 



^(^^e,^,fc _ f^sv'^^^'^ - v'^-^^^Sv'f^^ - 5v'^'"^^v'f^^)h^^ (s, y) 
xGe{t — s,x — y)dyds+ 

So /r3 (5*^^"^^ - 5'^'') (s, y)G^{t - s,x- y)dyds, 



(83) 
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where 

(• e,w,k . r I e,w,k \ / \ j \ r e,ui,fc 

/i,2't^,i + /i,3ti^,i + "f^i w,2j{y)dyj6v^^ 

X r e,w,k I f e,w,k \r \j \ e,w,k—l 

fi,20v{ +dv^^ w^2j{y)dyjvi^ 
/i,2'5t^^r''' + 5d'3"''''w;,2) (y)%)/i,3- 

(84) 

Next we observe preservation of polynomial decay of certain order in the 
schemes. This observation holds also for the naive iteration scheme v""'^, and 
it holds a fortiori for the iteration scheme v''''^'^ 

Lemma 5.1. Polynomial decay of order m of the functions v'^'^^'^ and 
their derivatives of order m is preserved if polynomial decay of 
order m is given for the functions v'^'^{t, .) G C"^nH"^, 5v'^'^{t,.) G C"^r\H"^ 
(for allt >0 ), and fi E C^^OH'^ and their derivatives of order m. An anal- 
ogous statement holds for the iteration scheme y'"'^'^^^ for the parameterized 
Cauchy problem. 

Proof. We prove the statement for the case m = 2 which is essential for our 
purposes. We can prove the lemma with the original source terms S^''^^^ 
for the original Cauchy problem which is is obtained by replacing w^2 by 
-^3,2(/i,i5 2^^) and wi by l2,,i{fi,i,v^2) Sti^^^ . The proof for the parme- 
terized Cauchy problem with the source terms Sw^^^ is analogous and sim- 
pler. The generalisation of the following argument to higher order m is also 
straightforward. We show that we have for m > 1 

\6v'^^+^\ < if > 1 (85) 

and similarly for We have to estimate convolutions with the Gaussian 

Ge. The expressions for 6v'^'''~^^ and ti'^-'^+i above are of the form 

/ h{y)G,{x - y)dy. (86) 

where h is some function which is a functional of v'''^{s,.) and 6v'''^{s,.) 
and fi are assumed to be in n C^, where the latter functions and their 
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derivatives up to esond order are assumed to be of polynomial decay of order 
m = 2. Furthermore we have the Gaussian factor G^it — s,x — .) is . We 
split up the integral of the convolution into two parts where one part is the 
integral for \y\ < On this domain we observe that the Gaussian has 
polynomial decay of any order. Indeed we have 



-s,x-y)\ = I 1=^ exp {-^) 



(2^e7r{t-s))" 

<\(t- s)™-"/2 1 { i^-v? X^ _A I eXD i 



(87) 



<\c{t-sr-/'j^l 



where 



^ _ I 1 1 f ix-y)^ \"' ( {x-yf 



is a constant depending on e, but finite for each e > 0. On the complemen- 

\x 
2 



tary domain \y\ > we need some properties of the integrand h. Next we 



observe 

c 



/{|,|>M}p?G,(x-2/)dy 
+ /||,|>M}&{|._,|>1,| - y)dy (89) 



for some generic constant C > 0. Hence the proof reduces to the observation 
that the integrand h in the form it has in the representation of Sv^^^'"^ is of 
polynomial decay. Since v""'^ and 5v'''^ are of polynomial decay inductively 
this implies that v""'^^^ is of polynomial decay too. Now upon inspection of 
we observe that 



-/i + + t;^'"'Vi,2)<5<f''^ + 5t;^'"''Vi,2t' ■ 

(90) 
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is of polynomial decay of order m = 2 by the inductive assumption of poly- 
nomial decay of order m = 2 of the functions v'^''^, 6v'''^, and v'^'''~^, and 
of their first order derivatives. The integrand (S"^'*^"*"^ — 5'^''^) (s, .) deserves 
more detailed consideration since it is itself a convolution. Note that for 

\x\ 

\y\ < we can use the Gaussian polynomial decay as above. Hence it is 
sufficient to estimate integrals of the form 



K4y-z)g{z)dz (91) 

{m>¥} 

where g is a functional of v''^'' ^v'^'''~^ and Sv'^''^ and their first order deriva- 
tives. Note that there are additional factors of the form u g'^ in the definition 
of ()84p but they have polynomial decay of order m by assumption and prod- 
ucts of functions of polynomial decay of order m certainly have polynomial 
decay of order m. Hence we may neglect this factor. Assuming 



< ^ for > T (92^ 



for some generic C > we can use a similar argument as above and write 
<\I[\y\>M]s,[\,\<M]KAy-z)9{z)dz\ 

+ 1 /{|,|>M}&{|,|>M} - z)9{z)dz\ ^93) 
+ 1 //u> M\^/i.i.M\ ^Ay - z)^dz\ go' ^ 



/||j,|>M}&{|,|>M} - z)^az\ t w yj^^. 

It remains to verify that g as implicitly defined in the expression for 5^^^ — S"^' 
satisfies the requirements needed for g in the last step. Now in ()84p all 

summands of the form ff^, (^^2^^ ) + ^ 2^^) t fi,2v''i are certainly of 

polynomial decay of order m as /i^2 is this by assumption and v'''^ and its 
derivatives up to order m = 2 by inductive assumption. 

Moreover, terms of the form fi,3l3,i{fi,i,v'^2^),v'^^^Li^2{fi,i,v%'') are of 
polynomial decay of order m since the integral forms are /3,i(/i,i, ^'^d 
are bounded l3,2{fi,i,v%^) and their factors /i^3 is of polynomial decay of 

order m by assumption and v""^ is of polynomial decay by inductive assump- 
tion. This shows that the considerations above apply for the first integral 
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in the expression for ()84p . Similarly for the other integrals in ()84p . We 
have shown that 6v'''^ is of polynomial decay of order m = 2 if v'^''' and 
dv"^'^ and their derivatives up to order m (and inductively all predecessors 
in the order of stage k) are of polynomial decay of order m. For the first 
order derivatives we have still locally integrable bounds of the first deriva- 
tives of the Gaussian such that the argument above can be repeated with 
standard a priori estimates for the first derivative of the Gaussian. For the 
second derivatives we shift over one derivative of the Laplacian kernel K^i to 
the integrand of form h or g above. These integrands involve derivatives of 
f,v''''',6v'''^ of first order. Hence the inductive assumption for order m = 2 
is still satisfied if these terms are differentiated. Furthermore we shift one 
derivative from the Gaussian to the Laplacian kernel K and then repeat 
the argument above concluding that second order derivatives of ti^-'^+i and 
j^<:,*:+i a,lso have polynomial decay of order m = 2. □ 

Next we prove a contraction property for the parameterized Cauchy 
problem which implies the existence of a classical solution v^''^ with v'''^{t, .) G 
n for all t > 0. Again w is a parameter in C|. 

Lemma 5.2. Assume that v'^'^'^ and y'^^^^f'-'^ and their derivatives up to 
second order are of polynomial decay of order m = 2. For any choice 

C > 4(2 + (2n + l)6Co)2C7i + 4((2n + 1)(C2 + (2n + 1)C3)(5C2 ^ 2C0C4) 
+2(^2 + (2n + l)C3)(2n + l)C7o8(2n + l)Co)) 

(94) 

we have 

jj ye,w,k ^e,w,k-i their derivatives up to order m are of polynomial 
decay of order m then a similar result can be obtained for an analogous norm 

Proof. We consider the case m = 2. The generalisation to m > 2 is straight- 
forward. We may rewite ()83p and write 

^(j^.,^,fc _ f^§y^^^^^ _ ^e,w,k^^^^,k _ g^e,w,k^^^,ky^ ,\^ y)dyds+ 
lo /iR3 (5'^'''^^ - Su,''^ (s, y)Ge{t - s,x- y)dyds 

_. ^yt,w,k+l,l _j_ ^yt,w,k+l,2 

(96) 
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and estimate the first integra term on the right side of ()96p separately, since 
it does not contain a Laplacian kernel, and can be treated differently there- 
fore. For the first and second order derivatives of Sv''''^''^'^^'^ we consider the 
representations 

xG^^i{t - s,x- y)dyds, 

(97) 

and 

^(^^.,^,fc _ f^Sv\^''' - v'''"'''5v'f''' - 6v'''"'^v'f''')fi,3) {s, y) 

) ,i 

xGe j{t — s,x — y)dyds. 

(98) 

For dimension n = 3 let 
C := 4(2 + {2n + l)6Co)2Ci + 4((2n + 1)(C2 + {2n + 1)C3)(5C2 + 2C0C4) 

+2(C2 + (2n + l)C3)(2n + l)Co8(2n + 1)Cq)) 

(99) 

(the reason with this definition will become apparent later). We assume 
inductively that 

In n \ 

sup \v^^^^\t,x)\+Y,\vT''^t,x)\ + J2 \vTi'(t^^)\ ^ Co, (100) 



i=l *=lii 



along with sup^g \v'^''^'^{t, .)\h'2 < C'o for the same generic constant Cq. In 
order to avoid inflation of the number of constants we assume w.l.o.g. the 
same upper bound for the relevant norms of the parameter function w, i.e., 

In n \ 

sup [\w{t,x)\+^\w,i{t,x)\+^\w,i,j{t,x)\\ <Cq, (101) 



t>0,x6M" 



1=1 i=l,i 



along with sup^>o \ w{t, .)\h^ < Cq. For the present we shall take this as in 
induction hypothesis at stage k (for A; = recall that we defined = — /i- 
hence, at stage k = Q the induction hypothesis is satisfied). Estimates in 
of Gaussians like Ge{t — s, .) or its first order derivatives can be obtained 
by breaking up integrals. Consider a function cpi £ C°° (i?i(0))), i.e. with 
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support in -Bi(O), and with (l)i{x) = 1 for |x| < 0.5. Note that (pi and 1 — (pi 
build a partition of unity on M^. We refer to [1], but this can be found 
in standard texts for sure. In addition we may use a generahzed Young 
inequality 

1/ < \f\Lp\g\Li, for ah 1 <p,q,r,< oo where i + i = 1 + i, 

(102) 

For i > s we have \Ge{t — s, < Ci < oo and |Ge,i(i — s, < Ci < oo 
for all 1 < i < n and some Ci > which is even independent of the size of 
t — s. We get the (generous) estimate 

exp{-Ct)\6v'''"'^+^'^{t, .)\h2 

< exp(-C7t) /o(2 + (2n + l)CofCi\6v'''"'^{s, .)\H2ds 

< exp(-Ct) /o exp(Cs)(2 + (2n + l)6CofCi\6v'''"'''f^]''^ds (103) 

^ (2+(2n+l)6Co)^Ci |^„,,.to.fc|exp,C7 
C 1"^ 1^2 

Hence, by choice of C > we have 

|^^.,«,,fc+i,i|exp,c ^ 1 .)| ^""2^'^. (104) 

Finally we do the estimate for i.e., the iJ^-estimate for 

/* / (S'/"-' - 5^') {s, y)G,{t -s,x- y)dyds (105) 
Using the generalized Young inequality this reduces to a if^-estimate for 

The latter term involves a convolution with the Laplacian kernel and we use 
the partition of unity in M.^ with (pi (support i?i(0)) and (1 — (pi) as defined 
above (cf. also [3] for similar considerations). For all 1 < i < n we have 



\(PiK,iiz)\d'z= / \(Piiz)\ YT^d-'z^ / \(Piiz)\ —r'dr = C (107) 
Jr3 7ir3 |r|^ 

for some constant C > 0. Hence, (piK^i G L^, and in order finish our 
prove of a contractive property we consider again three different parts of the 

exnrpssinn fnr Qk+1 _ ok _ ce,k+l,l _ ne,k,l ne,k+l,2 _ ne,k,2 nt.k+l,?, _ o<^,k,'i 
expiebbion iui Oyj — Ow n^Ow Ow n^Jw i^w 

considered above. All estimates can be done by splitting up the integrals 
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and considering the sums involving (j)iK^i G and (1 — 0i)J^ j € n L°° 
separately. We have finite constants 

C2 := \<PiK,,\li (108) 

and 

n 

C3 := 1(1 - </.i)i^,i|L2 + 1(1 - </)i)Ki|L- + ^ I ((1 - 4>i)K,i)j |l- (109) 

and we estimate the summands involving the (outside) truncated kernel 
(piK^i by the generalized Young inequality and the summands for the kernel 
(truncated inside) (1 — (pi)K^i with some basic techniques of the Fourier 
transform. In addition to the preceding inductive assumptions concerning 
the bounds for v'''^ and /i we mention our assumption on w. We have 

suPf>o.xeiR" '^{\wit,x)\ + I J2iW^i){t,x)\ 

(110) 

<C4 

for all t > and some constant C4. The first term 5'^'^'''"^'^ — Sw'''^ can be 
estimated by referring to the suprema. We remark that for bounded C^- 
functions g and functions 6v'''^ with bounded first order derivatives we have 
i?^-estimates for expressions of the form 

(111) 

for a constant C > which is an upper bound for \g\L°° and \gAL°°., i.e., 
all first order derivatives. In the following we fix time t > arbitrarily and 
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suppress the notation of time. This way we get 

fi,2v''^f''' + /i,3W^,i + f2,3W,2){y)dy]5v''f'''\^^ 



2 / , N 2 

<G 



,3 



III + (^^2"'''') + + ^^'2^) + /l,2^^'r''' + fl,3W,l + f2,3W,2^ 6v' 

+(2n + 1)C3| [fl + (y-f^'Y + + <f + /i,2<r'' + /l,3^,l + /2,3«^,2)5^^; 

< (2n + 1)(C2 + (2n + 1)C3)(5C72 + 2CoC4)|<57;;f ''^jj^i 

< (2n + 1)(C2 + (2n + 1)C3)(5C72 + 2C7oC4)|<57;^'"'''=|h2, 

(112) 

where we use the H^'^ bound for the factor (1 — 4))K^i and its first order 
derivative. For the second term we have 



3 fji 



I Qe,fc+1,2 _ Qe,fc,2| 



2 / , -,N 2 / , N 2 



(nx , e,w,k , e,w,k—l\ r- e,w,k . p c e,w,k , r t,w,k \/ \ t \ £,10,^—11 

+1 /r3(i - </'i)^,i)(a^ - y)( (^'^'^ + f'r'^"^) 

(n i" I e,w,k , e,w,k—l\ r e,w,k , j. r e,w,k , c- t,w,k \/ \ i \ 
2/1,1 + ^^,2 +"^,2 j '^^,2 +/i,2'5fj +<)t^,i w^2j{y)dyjvj^ 

< {C, + (2n + 1)C3)| ^3 ( {vT' + 

2/1,1 + ^^,2 +"^,2 j "^^,2 +/l,2')t',i +C't^,i ^^,2jjt^,i \h1 

< {C2 + (2n + l)C3)(2n + l)C7o8(2n + l)Co)|<5i;^'"''=(t, )\h^, 

(113) 

uniformly for ah i > 0, and where we suppressed the notation of the parame- 
ter t > for ah terms except the last one. The first term on the right side of 
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k-l\ 



the latter equation can be estimated using the generahzed Young inequahty. 

Note that the third term differs from the second term only with respect 
to the kernel function K2 instead of K^i. Hence, we have the same estimate 
for 

- ( /m3 K^2ix - y) ( {vf^ - (vf^ + + vf^ - + vf'^^ 

fl,2Svl'i + /l,3^3,l(/l,l! '^^2,2) + ^^,3 ^3,2(/l,l, f '2 ) + ^,3 " ^3,2 (/l, 1, ^^^^2 

im 

Summing up the estimates fr the source terms of the form Sw^~^^ — Sw'^ we 
have 

I 'Jw iJw I ^1 — 

< (2 + (2n + l)6CofCi + ((2n + 1)(C2 + (2n + 1)C3)(5C2 + ICod) 
+2(C2 + (2n + l)C73)(2n + l)Co8{2n + l)Co))\6v'''">'{t, .)\h2. 

(115) 

Hence, we observe that the choice 
C := 4(2 + (2n + l)6CofCi + 4((2n + 1)(C72 + (2n + 1)C3)(5C2 + 2CoCi) 

+2(C2 + (2n + l)C3)(2n + l)C7o8(2n + l)Co)) 

(116) 

is sufficient in order to get the |.|^^' with factor ^. □ 

We have obtained a solution v"^ of the parameterized Cauchy problem. 
Next a rather simple argument (using similar techniques as above) leads to 
a stability with respect to w. 

Lemma 5.3. The parametrized solution f is stable with respect to wi^t, .) G 
C| (uniformly w.r.t. t), i.e., for each e > we have 

.) _ y^,-^\t, .)\h2 < C\w{t, .) - W'{t, .)\c2 (117) 

for all t > and some constant C > 0. 

Now for each u; G we have obtained a solution f ^ of the parameterized 
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Cauchy problem which has the representation 

v^'-'it, x) = /o* /i,3 ((-/i + + v^'^h,2yf 

+iv''^-fiv:r'''-v^'^v%^)fi,3) 

+ /o ( kaz- y) {fi + (-r ) ' + (/i.i + <r) ' + 

- /o ( /m3 K2{x - y) {fl + (.r ) ' + + v:f) " + 

/l,2^^2 J + 11,3^,1 + v''fw^2^ {y)dyfl,3^ {S, z) 

X — . exp I — wT^\ I (izds. 

The divergence property of the vector field (/i, f^''^, w) leads to 

-W,3{t,x) = fl,l{t,x) 

+ /o7m3 ((-/l + /l/l,l+^^^'-/l,2Kr 

+ ^ye,^,k _ f^^e,n. _ ^^,-,yf )/i,3) G,,2(t -S,X- y)dyds 

+ /o ( /m3 ^,1 - y) {fli + (^r ) ' + + <f ) ' + 

/i,2w'i"' + /i,3t^,i + v'^fw,2^ ^)G'e,2(i - s, X - z)dzds 

- /o (/m3 ^,2(x - y){fl + (t^r)' + + <f )' + 

/l,2^^2,r + /l,3W^,l + ^,3"'w^,2) {y)dyf 1,3^ (S, z)Ge,2{t - S,X- z)dzds. 

(119) 

We can apply that f is stable with respect to w. A simple iteration scheme 
and contraction estimate leads to the divergence free solution (/i, u^''^, w) 
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of the Cauchy problem for each e > with w^^ G n H^. However we 
cannot ensure the L^-finiteness of the integral of w^^ with respect to the 
third variable over the whole space and therefore for the velocity V3 and /s. 
Applying this result we get 

00 

v^:=v^'' + Y.Sv^'''eHl^p^C- (120) 
k=l 

Hence for each t > we have 

v'{t,.) e CC" (121) 

for some a £ (0, 1) and uniformly with respect to time t. Hence we have the 
representation 

y'it,^) = /o /m3 {{-fl+flfl,l+v'fl,2)v% 

- fiv'f'" - v'v%)h,3)G,it -s-x- y)dyds 
+ !l ( ^3 K,^{z - y) (fl, + {v%)' + + v%)' + 

fi,2v\ + f 1,313 Afi^,^%) + f2,3h,2{fi,i^^%))iy)dy)v%){s,z) 

' ' ' (122) 

xGe(t — s;2; — z)dzds 

- /o ( /m3 K,2{x - y) [fl, + [v^^^f + (/i,i + v%f + 

h,2vl^i + /i,3^3,i(/i, 1,1^12) + v%l3,2{h^^v%)){y)dyh;i){s,z) 

xGe{t — s; X — z)dzds 

of the solution v'' of the Cauchy problem ()68p . Spatial regularity can be ob- 
tained from this expression straightforwardly then. Furthermore, the argu- 
ment above can be repeated for any order m > 2 of regularity /i G n 
yielding v''{t, .),w{t^ .) G C^CiH^ uniformly with respect to t and with poly- 
nomial decay of order m. Furthermore, using standard a priori estimates 
for the Gaussian as in the lemma on polynomial decay above we observe 
from (|122p we see that \v'^{t, .)\f{2 < C for some constant C > which is 
independent of e. Furthermore you easily observe that this holds also for 
the solution := limji-oo v'^{t, .) of the related elliptic problem. 

^kampenOmathalgorithm . de , kampenSwias-berlin . de. 
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^In this paper we did not state explicitly that we can derive linear growth in time of the 
squared value function integrated over space (and linear growth with respect to time of 
the derivatives up to second order of the squared value function integrated over space). We 
said that we have linear growth of v(/, .) with respect to the time step number / and with 
respect to the _H"^-norm maxi<i<„ \vi{l, ■)\h^ is sufficient in order to get a global scheme 
although the Leray projection terms contains squared (first order derivatives of integrals 
of components of the) value function. However, it is easy to show that the scheme leads to 
linear growth with respect to the time step number and the squared T^^-norm \vi{l, .)|^2 
for the components of the value functions, and hence is an alternative to the controlled 
systems in [1] and [2]. In a third part of this paper we shall discuss such aspects in more 
detail. 
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